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1 Introduction 

Quasilinear hyperbolic systems have a special place in the theory of partial dif- 
ferential equations since most of the PDEs arising in continuum physics are 
of this form. Well-known examples are the Euler equations for a perfect com- 
pressible fluid, the equations of elastodynamics for a perfect elastic solid, and 
equations describing a variety of field-matter interactions, such as magnetohy- 
drodynamics etc. It is well-known that for all these systems the Cauchy problem 
is well-posed, i.e., it has a unique classical solution in a small neighborhood (in 
space-time) of the hypersurface on which the initial data are given. 

On the other hand, it is not expected that these systems will have a global- 
in-time regular solution, because shock discontinuities are expected to form at 
some point, at least as long as the initial data are not very small. In more 
than one space dimension, there are no general theorems to that effect however, 
mainly because in higher dimensions, the method of characteristics, which is a 
powerful tool in one dimension for the study of hyperbolic systems, becomes 
intractable. 

In 1985 T. C. Sideris published a remarkable paper on the formation of 
singularities in three-dimensional compressible fluids , proving that the clas- 
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sical solution to Euler equations has to break down in finite time. His proof was 
based on studying certain averaged quantities formed out of the solution, show- 
ing that they satisfy differential inequalities whose solutions have finite life-span. 
Such a technique was already employed by Glassey || in the case of a nonlin- 
ear Schrodinger equation. The idea is that by using averaged quantities one is 
able to avoid local analysis of solutions. The same technique was subsequently 
used to prove other formation of singularity theorems: for a compressible fluid 
body surrounded by vacuum in the nonrelativistic @ and relativistic |l(J cases, 
for the spherically symmetric Euler-Poisson equations in the attractive and 
repulsive ]j| cases, for magnetohydrodynamics ||], and for elastodynamics Q. 
In this paper we present two more such "Siderian" blowup theorems: one in 
relativistic fluid mechanics, and the other in plasma dynamics. 



2 Relativistic Fluid Dynamics 

Let (M,g) be the Minkowski spacetime, with (x' 1 ), p = 0, . . . , 3 the global 
coordinate system on M in which g^ = diag(— 1, 1, 1, 1). We will use the 
standard convention that Greek indices run from to 3, while Latin ones run 
from 1 to 3. Indices are raised and lowered using the metric tensor <?, and all up- 
and-down repeated indices are summed over the range. We also denote t = x a 
and x = (x 1 , x 2 , x 3 ). In the following, we adopt the notation and terminology 
of [Q and quote from it some of the basic facts regarding relativistic dynamics: 
The energy tensor for a relativistic perfect fluid is 

= {p + p)u^u" +pg^. (2.1) 

In this formula, 

1. u = (u°,u) is the four- velocity field of the fluid, a unit future-directed 
timelike vectorfield on M, so that g(u, u) — —1 and hence 

u° = Vi + H 2 . 

We note that here, unlike the nonrelativistic case considered by Sideris, 
all components of the energy tensor are quadratic in the velocity. 

2. p > is the proper energy density of the fluid, the eigenvalue of T cor- 
responding to the eigenvector u. It is a function of the (nonnegative) 
thermodynamic variables n, the number density and s, the entropy per 
particle. The particular dependence of p on these variables is given by the 
equation of state 

p = p(n,s). (2.2) 

3. p > is the fluid pressure, defined by 

P = n%- p. (2.3) 
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Basic assumptions on the equation of state of a perfect fluid are 

^ > 0, ^ > 0, ^ > and = iff s = 0. (2.4) 
on on Os 

In particular, these insure that n, the speed of sound in the fluid, is always 
real: 

2 ( d V S 



\dp 

In addition, the energy tensor ( [2.1| ) must satisfy the positivity condition, 
which implies that we must have 

p<p. (2.5) 

A typical example of an equation of state is that of a polytropic gas. A perfect 
fluid is called polytropic if the equation of state is of the form 

p = n + S?L n i (2.6) 
7-1 

where 1 < 7 < 2 and A is a positive increasing function of s (The speed of light 
is equal to one). This implies that p = An' and thus the sound speed n(n, s) is 
determined as follows: 

dp\ dp/dn 7(7 — l)^4n 7_1 



V = -7- = 



dp J dp/dn 7 — 1 + 7/ln 7 1 
In particular, the sound speed is increasing with density and is bounded above 

by V7- !• 

The equations of motion for a relativistic perfect fluid are: 

d v T^ = 0. (2.7) 
Moreover, n = n(x) satisfies the continuity equation 

d„(nu") = 0. (2.8) 



Given an equation of state ( |2.2[ ), the system of equations (2.7-2.8) provides 5 
equations for the 5 unknowns n{x), s(x) and u(x). The component of (2.7) in 
the direction of u is 

u"d v p+ {p + p)d u u u = 0. (2.9) 
As long as the solution is C , this is equivalent to the adiabatic condition 

u v d v s = 0. (2.10) 



The component of (2.7) in the direction orthogonal to u is 

{p + p)u v d v u fi + h> J - v d u p = Q, (2.11) 
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where 

is the projection tensor onto the orthogonal complement of u(x) in T X M . 

Thus the system of equations for a relativistic fluid can be written as follows: 

d v (nu v ) = 0, 
{p + p)u v d v vr + bP u d v p = 0, (2.12) 
u v d v s = 0. 

The Cauchy problem for a relativistic fluid consists of specifying the values of 
n, s and u on a spacelike hypersurface Eo of M , 

n\ Tjn = n , si„ = s , ui En = u , (2.13) 



and finding a solution (n, u, s) to (2.12,2.13) in a neighborhood of Eo in M. In 
particular, let Eo = R 3 x {0} be the hyperplane t — in M and suppose the 
initial data (2.13) correspond to a smooth compactly supported perturbation of 
a quiet fluid filling the space, i.e., assume 

no, so and Uo are smooth functions on R 3 and there are positive 
constants Ro, n and s such that outside the ball Br o (0) we have (2-14) 
no = n, so = s, and Uo = 0. 

Let fj — rj(n, s) be the sound speed in the background quiet state. We then have 



PROPOSITION 1. Any C 1 solution of (fDj|Dj,(DI|) will satisfy 

n — n, s — s, u = 0, 
outside the ball B R ^(Q) where R(t) = Rq + fjt. 



Proof. It is enough to check that the system (2.12) can be written in symmetric 
hyperbolic form: 

A^{U)df,U j = where = Af t and A% is positive definite. (2.15) 

This can be accomplished for example by using p instead of n as an unknown. 
By (2.4), we can think of n as a function of p and and s and thus of p as a 
function of p and s. By (2.3) it is then easy to see that (2. IS) is equivalent to 
the following system for the unknowns U = (p, u, s): 



-u v d v p + nd v u v = 

KP+PM 

nh^dvp + (p + p)nu v d v u^ = 

u v d v s = 0. 



(2.16) 
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Let U — (p, 1, 0, 0, 0, s) denote the constant background solution to (2.16). Let 
C := (p + p)fj > 0. We then have that the differential operator P = A^d^ 
corresponding to the linearization of (2.15|) at U is symmetric hyperbolic, with 



A°(U) = diag( , , C, C, C, C, 1) , A 1 = A\U) = 



/ fjef \ 



fjei 

\ / 



Once we have this, we can use energy estimates, as in Jjjj, to conclude the 
desired domain of dependence statement. 



We now prove that for large enough initial data, the solution to ( 2.12 , 2.13 2.14 ) 
cannot remain C 1 for all t > 0. Such a result was announced in |lC(|, but the 
unpublished proof contained an error which invalidated the argument |p"l| . 

First of all, a scaling analysis shows that without loss of generality we can 
set Ro = 1. Let 

B t = {x e M | x° = t, |x| < R(t) = l + fjt} 
denote the time t slice of the range of influence of the data, and let 



Q(t) 



9ij' 



irpOj 



J x • uu°(p +p) 



(2.17) 



be the total radial momentum of the fluid at time t. We then have 

Q'{t) = j I j ihi 1 ' 



Let 



E 



T 



on 



r 



00 



(p + p)\u\ 2 + 3(p-p). 
{p + p)\u\ 2 + p - p 



(2.18) 



(2.19) 



be the total energy of the perturbation. By (2.7) it is a conserved quantity, 
E = Eq. Our goal is to use E to obtain a differential inequality for Q that 
would lead to blowup. 

We are going to make two assumptions on the equation of state of the fluid, 
which are quite natural from a physical point of view. First we note that, as 
mentioned before, we can use the pressure p as a thermodynamic variable in 
place of n. The equation of state of the fluid then has the form p = p(p, s). The 
two assumptions are: 

(Al) p(p, s) is a non-increasing function of s, for each p. 
(A2) 7y(p, s) is a non-decreasing function of p, for each s. 
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These two assumptions are in particular satisfied for a polytropic equation of 
state (2.6). In order to see that, we observe that n = [p/A{s)) 1 ' 1 , and from 
there we get 

1 1 



l/ 7 



It is then clear that (Al) holds. Moreover 

7(7 - l)AV7( s ) p (7-i)/7 



7 - 1 + 7A 1 /7( S )j)(7-l)/7 



shows that (A2) is satisfied. 

We also make the following assumptions on the initial data: 

(di) n<\- 

(D2) E > 0. 

(D3) s (x) > s for all x e B Q . 

By ( |2.10D , the entropy s is constant along the flow lines, and thus (D3) implies 
that s(x) > s for x £ Bt- By (Al) and (A2) we then have 

P-P = P(P, s ) - PiP, s ) = P(P, s ) - PiPi + p(p, s) - p(p, s) < pip, s) - pip, s) 



-~-ip ,s)dp 
5 dp 



p rfip',s) dp rj 2 ^ 



By ( |2.18| ) and ( |2.19| ) we then obtain 



Q'it)>ifE+{l~if) Kp + p)W\ 



which implies, by virtue of (Dl) and (D2) that 

Q'(t) > (1-377 2 ) J ( p + p)\u\ 2 >o. 

In particular Q{t) > if Q(0) > 0. 



On the other hand, we can always estimate Q(t) from above, using (2.5): 



Q\t) < ( /( P + 7J )|u| 2 )i? 2 (t)( / (p + p)(|u| 2 + l) 

B t 



< 2 [ {p + p)\n\'\R\t)[ (p + p)\u\'+p-p + p 



< 



-2-^Q<(t)I?{t)[E+^pR\t)\. 



Integrating this differential inequality and changing the integration variable 
to r = R{t), we obtain 



1 1 1 - 3r7 2 



2 r R(t) 



Q{t) ~ Q(0) 2r7 



dr 



Er 2 + ^-pr 5 ' 
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which contradicts the positivity of Q for all time provided the initial data sat- 
isfies the following final assumption: 



(D4) Q(0) > 



2fj 



1 - 3jf 



dr 



Er 2 



±§-pr 5 



The contradiction implies that there exists a certain T* < oo by which time a 
C 1 solution has to have broken down. In particular, the domain of dependence 
may break down at an earlier time, perhaps because a shock discontinuity forms. 
We have thus proved 

THEOREM 2.1. Suppose that the equation of state of a fluid satisfies (Al) 
and (A2). Then the Cauchy problem (2.12, pTi~3 2.14) with initial data satisfying 
(Dl— D4) cannot have a global-in-time C 1 solution. 



Remark 1. It is easy to obtain a simpler, sufficient condition for blowup: Let 

f(y) ■■= 



dr ' 1 



r 2 (r 3 + y) 



By (D4) we thus need 



It is easy to see that /(0) = 4, /'(0) = 16/7 and that / is a concave function of 
y, so that f{y) < y-y + 4- It is therefore enough to have 



Q(°)> ^ ?)■ ( 2 - 21 ) 



7(1 -zn 2 y 3 

We note that unlike the nonrelativistic case, the lower bound for the initial radial 
momentum in (D4) or (2.21) depends on the initial energy, and thus on the 
initial velocity. Since Q is of the same order of magnitude as E, it is worthwhile 
t o sho w that there exist data sets satisfying these largeness conditions. In fact, 
( [2.21 ) can be satisfied for n small enough. All that is needed is dp/dn > at 
n — 0. We illustrate this in the following by considering the poly tropic case. 

Let us consider a fluid with a polytropic equation of state (|2.6|) , and consider 
initial data of the following form 

no(x) = ni/>(r), u (x) = -<f>(r), s (x) = s + 0(r), (2.22) 

r 

where r = |x|. 4> an d ip are smooth, positive functions on [0, oo) such that 

4>(r) = for r > 1, ^(0) = 0, (2.23) 

and 

ip(r) = 1 for r > 1, / (i/>(r) - l)r 2 dr = 0. (2.24) 
Jo 
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We then compute 



E 



(po +Po)|u | 2 + po - p 



Ann 



7-1 



57-1 



[A{s)^(n(t) 2 + 1) - A(s)} \ r 2 dr, 



and thus E > by (D3) and Q2.24| ). Now, 



Q(Q)=47rn/ <h/l + <W + ^ 



(i 



7-1 



Dividing (2.21) by 47rn, all we need is that the following inequality be satisfied 
for n small enough: 



32 



' VWl + </> 2 r 3 dr+0(n 7 " 1 ) > — Tr^fjl ^4> 2 r 2 dr + — + O^" 1 ) 

7(1 -3?7 2 ) t7 12 

(2.25) 

This is clearly true since fj — > as n — > 0. We have thus shown 

PROPOSITION 2. Lei and "0 oe (wo smooth, positive functions on [0,oo) 
satisfying ( 2.23| ,2.24). Then there exists n > small enough (depending on <f>, 
tp and 7 J such that the initial data set (no, Uo, sq) of the form (|2.22| ) satisfy the 
conditions (Dl— D4), and thus lead to a blowup for (2.12 ,pTl3[ 2 . 1 < 



3 Euler-Maxwell with Constant Background Charge 

A simple two-fluid model to describe plasma dynamics is the so called Euler- 
Maxwell system, where a compressible electron fluid interacts with a constant 
ion background. Let n(i,x), s(i,x) and v(t,x) be the average electron density, 
entropy, and velocity, let n be the constant ion density, and let E(t, x) and 
B(i,x) be the electric and magnetic fields. Let c = speed of light in vacuum, 
e = the charge of an electron, and m = the mass of an electron. The Euler- 
Maxwell system (see || pp. 490-491]) then takes the form: 

d t n + di(nv l ) = 

W+diT> - —E ^._ c(VxB ). = -Anenv\ (3 ' 1} 

d t s + v l diS = 

together with the constraint equations 

d i E i = Ane{n - n), d i B i = 0. (3.2) 

In the above, II is the momentum vector, 

n = nvH - — (E x B), 

Annie 
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and T is the stress tensor, which can be decomposed into material and electro- 
magnetic parts: T = Tm + Tk, with 



Til = rwV + — P S ij , 
m 



12? = — f-(|Er + mns 13 -e'e 1 -b 1 b 3 \. 

p is the electron pressure, which is modeled by a polytropic law p(n, s) — A{s)n 1 , 
where 7 > 1 and A is a positive increasing function. 



The system (3.1) being hyperbolic, we once again have the domain of de- 
pendence property. However, this time the largest characteristic speed in the 
background will be c, the speed of light. We recall that in Sideris's original ar- 
gument , the largeness condition on the initial data implied that the initial 
velocity had to be supersonic at some point, relative to the sound speed in the 
background. An analogous result in the Euler-Maxwell case would thus require 
that the initial velocity be superluminar at some point, which is absurd. How- 
ever, we note that if the data is spherically symmetric, so will be the solution, 
and thus there will be no electromagnetic waves, and the largest characteris- 
tic speed will once again be the sound speed, so a Siderian blowup theorem is 
possible in the spherically symmetric case. Moreover, since spherical symmetry 
implies that the flow is irrotational, such a blowup result is complementary to 
the recent construction |i) of global smooth irrotational solutions with small 
amplitude for the above system. We note that a blowup result in the spheri- 
cally symmetric, isentropic case with no background charge has been obtained 
[Pj using Riemann invariants. 

Remark 2. Under the assumption of spherical symmetry, the Euler-Maxwell 
system reduces to what is often referred to as the spherically symmetric Euler- 
Poisson system (with repulsive force). We note the important distinction be- 
tween this, and the general Euler-Poisson system obtained by taking the New- 



tonian limit c — > 00 in (3.1). The latter is not a hyperbolic system, and does 



not have finite propagation speeds. 

We have the following theorem: 
THEOREM 3.1. Let vq, <7o and Uq be smooth functions on R + satisfying 

u {r) = a Q {r) = v (r) = for r > 1, u (0) = 0, <r (r) > 0, 
and the neutrality condition 



1 

2 



v (r)r z dr = 0. (3.3) 



i) 



Let s > be fixed. Then 
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(a) There exists T > and functions v, <r,u, E £ C 1 ([0, T) 



such that 



i/(0, r) = v (r) 
<r(0,r) = <7 (r) 
u(0,r) = uo(r) 



Aire 



vo{r')r' dr' . 



and such that the Euler- Maxwell system (3.1) has a unique solution of the 
form: 



n(t, x) 
s(t,x) 



n + u(t, r) 
s + o~(t, r) 



v(t,x) = u(t,r)~ 
r 



E(t,x) = E(t,r)- 
B(t,x) = 0, 



(3.4) 



where r = |x|. 

(b) Fort G [0, T), (n, s,v, E) satisfy the reduced Euler-Maxwell system: 



d t n + di{nv l ) 
d t s + v l diS 

dtinv') + djT ij 
d t E l + Airenv 1 






en . 

m 
0. 



-J5* 



(3.5) 



where 



rptj 



1 



1 -(-\E\ 2 S z: > - E l E j ), 



together with the constraint Poisson equation: 

d t E l = 47re(n-n). 

(c) Let fj — \J 7^4(s)nT _1 6e i/ie sound speed in the background, R(t) 
and let 

Dr :={(*, x) | < t < T, |x| > R(t)}. 



\ + fjt, 



Then we have (n, s, v, E) = (n, s, 0, 0) on Dt- 

(d) For any fixed vo(r) which satisfies ( |3.3[ ) ; there ex ists Uo(r) sufficiently 
large, such that the life-span of the C solution ( |3 .4[) is finite. 



Proof. (a) The Euler-Maxwell system ( |3.1|) can be written as a positive, sym- 
metric hyperbolic system, and therefore has a unique, local C 1 solution 
with n > provided its initial data are sufficiently smooth. Notice that 
the initial data are spherically symmetric. Because of the rotational co- 
variance properties of the Euler-Maxwell system and the uniqueness of the 
local solution, the solution remains spherically symmetric and (a) follows. 



(b) follows since B i = 0. 
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(c) Notice that from the Poisson equation at t = 0, 

47re 



£7(0, r) 



vo(r)r dr = for r > 1 



by th e ne utrality assumption (3.3). Now the reduced Euler-Maxwell sys- 
tem (3.5) is still a hyperbolic system, and we can deduce (c) via the 
Proposition in p2[. 



(d) Let 



Q(t) := ±- / a: • n = / 



A direct computation yields: 



Q'(t) = 



rax + —(p-p) + 

TO 07TTO 



rnu r dr. 



2 dr 



r£7 r 2 dr, 



where p = p(n, s). Meanwhile, by the first and fourth equations in (ETa), 



rE(t,r)r 2 dr = / r£7(0, r)r 2 dr - Ave / Q(t')dt' 
i) Jo Jo 



Integrating by parts, we notice that 



r£7(0, r)r dr = — 47re / ^o(r)r <ir. 
(i Jo 



We now define j/(t) := / Q(t')dt' and obtain 
Jo 

y"(i)+a, 2 y(*) = G(i), 

2 Aite 2 n . 



(3.6) 



where or 



TO 



is the plasma frequency, and 



G(t):=-Lu 2 v Q (r)r i dr+ I nu H (p — p) + —^—E 2 \ r 2 dr. 

' Jq \ TO 87TTO 



Therefore, form solving the ODE (3.6) for y(t), we have 



y"(t) = -uy'(0) smut + G(t) - oj sinu{t - T)G(r)dT 



(3.7) 



We recall the conserved quantities energy: 



1 



1 2 

— nu H ; r 

2 to( 7 -1) 



(A(s)n 7 - A(s)fi 7 ) + £7 2 ^ r 2 dr, 

07TTO 



and mass 



M = — / (n-n) = v(t, r)r 2 dr. 
4?r Jhs Jo 
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From the neutrality condition ( |3.3| ) we have M = 0. Also, s(0, x) > s since 
(To > by assumption. By the adiabatic condition (the second equation in 
(3.5)) entropy is constant along flow lines, and thus s(i, x) > 5 for t <T. 
Since A is an increasing function, 

A(s)n 7 - A{s)n« > A{s) J rf - ff > f} 2 J n-n = Q. 

Hence we have 

/>oo 

a£ < G(t) +lu 2 v (r)r 4 dr < {3£, 
Jo 

with a = min{l,3(7 — 1)}, j3 = max{2,3(7 — 1)}. But for large enough 
u (r), J °° v (r)r 4 dr is dominated by £(0). Hence, we have 

|£ < G(i) < 2/3£ 

for sufficiently large uo(r). Moreover, we have 



Q (t) < R (t) / nu 



n < CR b {t)n£. 



(3.8) 



C will henceforth denote a generic numerical constant. By choosing Uq(t) 
large such that £ (t) = £(0) > 1, we have 

|y'(0)| = |Q(0)| < cVn£. 



Thus from (^J), there exists T = T (7, n,u)>0 such that for < t < T , 

Q'(t) > C£. 



Together with (3.8), we deduce for < t < T , 



Q'it) > 

Integrating over [0, T ] we obtain 



C 



R 5 (ty, 



-Q 2 (t). 



1 



Q(0) Q(T )-nfj l (1 + fjTo) 4 



(3.9) 



We can the n c hoose uo(r) sufficiently large, so that Q(Q) is so large to 
contradict (3.9). 
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